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Abstract 

For the first time, a phenomenological analysis of the experimental electromag- 
netic form factors of the nucleon, both in the timelike and spacelike regions, is 
performed by taking into account the effects of nonvalence components in the nu- 
cleon state, within a light-front framework. Our model, based on suitable Ansatzes 
for the nucleon Bethe-Salpeter amplitude and a microscopic version of the well- 
known Vector Meson Dominance model, has only four adjusted parameters (deter- 
mined by the spacelike data with x^/ datum ~ 1.7), and yields a nice description of 
the experimental electromagnetic form factors in the physical region in the range 
—30 (GeV/c)'^ < < 20 (GeV/c)'^, except for the neutron one in the timelike re- 
gion. Valuable information can be gained in the timelike region on possible missing 
Vector Mesons around ~ 4.5 (GeV/c)^ and q'^ ~ 8.0 (GeV/c)'^. 
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In recent years there has been a renewed interest in the investigation of nu- 
cleon electromagnetic form factors (FF), given an unexpected discrepancy 
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between experimental data for the spacelike (SL) ratio /i^ G^^{q^) / G\^{q^) ex- 
tracted by using: i) the Rosenbluth separation method (see Ref. [1] for recent 
measurements) and ii) the polarization transfer technique adopted in experi- 
ments carried out at TJLAB [2]. Indeed in the SL region (where the squared 
four- momentum transfer becomes negative, i.e. = — |qp < 0), data ob- 
tained by the Rosenbluth separation follow the dipole law, while, surprisingly, 
data from the polarization transfer technique decrease faster than the dipole 
law for = —q^ > 1 (GeV/c)^. This experimental puzzle has not yet been 
completely explained, although both two-photon exchange processes [5] and 
higher-order radiative corrections [1] appear relevant for its solution. 

Furthermore, the timelike (TL) region calls for both experimental and the- 
oretical investigations (in particular for the neutron !), since the ratio, R, 
between experimental neutron and proton form factors, beyond the thresh- 
old = 4M^ (with the nucleon mass), turns out to be greater than 
one [3], while naive expectations from perturbative QCD (see, e.g. [6j) yield 
R ~ \ed/eu\ = 0.5. A deeper understanding of all these experimental issues 
could open new windows in the investigation of the nucleon internal structure, 
also elucidating the role of small components in the nucleon state, (see, e.g., 
[7j for the possible influence of the above mentioned SL puzzle on the nucleon 
shape studies). 

Within the light-front dynamics [8.9J, we successfully reproduced the pion 
experimental FF in the interval -10 {GeV/cy < < 10 {GeV/cf [TO] , 
namely both in the SL and TL regions, by introducing components of the 
pion state beyond the valence one. Aim of this letter is the generalization of 
our approach to the nucleon, presenting for the first time within the light-front 
dynamics a unified, direct calculation of both SL and TL nucleon FF (see, also 
jll] for preliminary results with a slightly different model). In particular, the 
role of the contribution due to the gg-pair, created by the incoming virtual 
photon, turns out to be essential as in the pion case. Our approach shares many 
ingredients with the model of Ref. [12], but it exhibits distinct features, like 
i) the choice of a reference frame with the plus component of the momentum 
transfer g"^ = g° + gz 7^ [13], allowing a unified analysis of SL and TL 
regions (see Figs. [1] and [2] for a diagrammatic illustration), and ii) the gauge- 
invariant dressing of the quark-photon vertex through a microscopic Vector 
Meson Model (VMD) flO]. 

Nucleon FF, that enter in the macroscopic description of the em current op- 
erator, I^{q'^), are calculated in a reference frame where = Pn± = and 
g"*" = \J\q'^\- In the SL region, where q^ = — (with P/v and Pj^ the 
initial and final nucleon four-momenta, respectively) and g^ < 0, the nucleon 
Sachs FF are given by 
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where rj = — 2M^/g^. The expressions for the TL form factors corresponding 
to Eq. ([T]) can be easily obtained by changing P^r with —P^ and with P/y. 
In our approach, for the SL kinematics, the matrix elements of the nucleon 
current (see [[11]) are approximated microscopically in impulse approximation 
by the Mandelstam formula [H] as follows (see also [12]) 



S-\h)I^{h,q) <l>'^{h,h,h,PN)} (2) 

where is the nucleon Dirac spinor, the factor 3 comes from the symmetry 
of our problem, is the number of colors, {ki,k2-,k^i Pn) the nucleon 
Bethe-Salpeter amplitude (BSA), fcj the i-th constituent quark (CQ) four- 
momentum, k'^ = k^ + q, P/v = ki + k2 + ks and = ki + k2 + k'^ . In Eq. ([2]), 
J2 implies a sum over isospin and spinor indexes, S{ki) is the Dirac propagator 
of a CQ (with a chosen constituent mass m = 200 MeV) and X^{k^, q) is the 
quark-photon vertex, obtained by dressing a pointlike quark (see below). An 
expression analogous to Eq. ([2]) holds for the TL region (see [T5]). 

The nucleon BSA must have a Dirac structure, that has been devised exploit- 
ing a qqq — N effective Lagrangian which couples a scalar-isoscalar quark-pair 
plus a quark to the nucleon, as suggested in Ref. [12j. In particular, for the 
present calculation, no derivative coupling has been considered. Then, the 
properly symmetrized BSA of the nucleon is approximated as follows 

$^^(^1, k2, ks, Pn) = A(A;i, k2, ks) W{ki, k2, ks, Pn) (3) 

where A{ki, /c2, k^) describes the symmetric momentum dependence of the ver- 
tex function upon the quark momenta, and W^is given by 

= [5(123) + 5(312) + 5(321)] Xr^ U'^n (4) 
with Xtat the nucleon isospin state and 



S{h) Ty 7^ Scik,)C ® Sih) 



(5) 



In Eq. ([5]), C is the charge conjugation operator, Sc{k) = C S^{k) C ^ and 
the symbol ® keeps separated the matrices acting on the quark pair from 
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the ones acting on the quark-nucleon system. Note that the symmetry in the 
quark pair variables reduces the number of possible terms from 6 to 3. 

The quark-photon vertex X^{k,q) = Xjg + Xjy has both isoscalar and 
isovector contributions. In turn, each term, Xf (with i = IS, IV) , is the sum 
of i) a purely valence (V) contribution, I-'y{k,q), which is present in the SL 
region only (see below) and ii) a nonvalence (NV) contribution, Xf^y(/c,g), 
corresponding to the gg-pair production (Z-diagram). Moreover, I^i^y{k,q) is 
composed of a point-like bare term and of a VMD term (cf [16]). Summarizing 
one has 



^NV% q) = e{q+ + k+) e{-k+) [Zb M 7^ + Z{.^,r^^^(fc, g, z)] (6) 

where Mis = 1/6, AZ/y = 1/2. The constants Zb, Zy^j and Zy\j are unknown 
weights for the pair-production contributions, to be determined from the phe- 
nomenological analysis of the experimental data. In principle, we should ex- 
pect only one renormalization factor, but in the actual fitting procedure, we 
have taken Zy\i ^ ^v"m = ^b, given the lower degree of knowledge of VM 
isoscalar sector. We anticipate that, from our fitting procedure, the devia- 
tion from the equality is ~ 10%. As in the case of the pion [10], the bare 
term fulfills the current conservation in a covariant model [T5]. For the 
VMD term, we extended the microscopic model of Ref . [10] , by including the 
isoscalar mesons and by making the VM vertex trivially transverse to q^ [15] 
(this means q-TyMD = 0). The same VM mass spectrum, em decay constants 
and total decay widths of Ref. [TU] have been used for the isovector part of 
the VMD term. As to the isoscalar term, for i = 1,2,3, VM masses and the 
corresponding total decay widths have been taken from PDG [17J, while for 
i > 3 we have calculated the masses by using the mass operator of Ref. [18], 
with an interaction parameter wjs = wjv — 0.27 GeV^ {wiv = 1-556 GeV'^), 
in order to follow the Anisovitch-Iachello law (see, e.g., [TU] for the isovector 
case), and we have adopted the same total decay width = .150 GeV as we 
had for the isovector case. The em decay constants, F*+g_, necessary for de- 
termining Fy^^(fc, q, i), have been calculated with the model of Ref. [18], and 
agree, within the errors, with the corresponding experimental values of the 
known IV and IS vector mesons [T7]. Finally, we considered up to 20 mesons 
for achieving convergence at high 

Following the pion case [TU], the four-dimensional integrations on ki and k2 
in Eq. are regularized by assuming a suitable fall-off of the momentum 
component of the BSA. Furthermore, in the integrations on fcf and A;^ we 
consider only the poles of S{ki), namely we disregard the analytic structure of 
A(A;i, k2, ks) and of the momentum components of the VM amplitudes, present 
in Tyj^j-,{k, q, i) (see [ID]), since it affects Fock sectors beyond the ones implicit 
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in Figs. [T] and [2l Then the covariance is only approximate (see Ref. [TO] for a 
quantitative discussion in the pion case). 

For the sake of concreteness, let us show the formal expression of the micro- 
scopic current I^{q^) (whose matrix elements are given in Eq. ([2])) in the SL 
region. It becomes the sum of two contributions: i) a purely valence (or tri- 
angle) contribution, Iy[SL,q'^) (Fig. [H diagram (a)), where both the nucleon 
vertexes have two quarks on their A;~-shell {k~ = = (|k_|_p + m^)//c+) and 
the quark variables are in the valence region [P^ > kf > 0); ii) a nonvalence 
(pair-production or Z-diagram) contribution, I'^y{SL,q'^), where the initial 
nucleon vertex has a quark outside the valence range {k^ < 0, see Fig. [H 
diagram (b)) viz 



I^{SL,q' 



3 Nr 



dkt 

2(2vr)6 J 14 



dk2 f dki^ d]i.2± 

U+ 1,1+ 



^<*^{k^M,P'N) X 

(7) 



I'NviSL.q^ 



3 



dkt 



2(27r)6 7 kt J ^ k2 



dkt f d ki± dk2± ^^j^^ 
kt k'+ 



e{kt) X 



^3 ^^3 



9 ^'son + ^3on 



(8) 



where ki = {kf, kj_|_} is the light-front momentum and k'^^^^ = {P'j^ — ki — k2) 
The quantity JF'^ is a 4 x 4 matrix (see [12] and [15]) constructed from and 
Xf(A;,g) (Eq.ED, given by 



JTM = (1^^ + ^) jM (1^3 + ^) Tr /C(2, 1) 

+ {^'^ + m) Xj^(^3 + m)/C(2, 1) + /C(2, + m) Xj^(^3 + m) 
-/C(l, 3) 7'[X,^5 - X^y]7'(^3 + + m) 

+2(^2 + m) Tr0^ + m)Jj^5/C(3, 1)] (9) 

where = 2js + tn 1'iv ^-nd }C{i,j) = + m){f,j + m). 

In Eqs. dZl E]) the momentum dependence of the vertex functions in the va- 
lence range (P^ — — 0) is expressed through a light-front wave function, 
\&Ar, which is a PQCD inspired wave function a la Brodsky-Lepage (see, e.g., 
[8]), described in terms of the squared free mass of the three-quark system 
M^f^{ki, k2, ks) = P^ J2i Kon: i-e- 
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N 



2, , 



(9 m 



2\7/2 



7/2 



(10) 



where = kf / and A/" is a normalization constant, obtained from the 
plus component of the proton current at = 0, i.e. from the proton charge 
normalization. In Eq. fIlOp the power 7/2 and the parameter p = 0.13, which 
controls the end-point behavior and affects the FF mainly through the Z- 
diagram, allow one to obtain an asymptotic decrease of the valence contri- 
bution faster than the dipole G^dg^l) = [1 + \q'^\/i0.71{GeV/cf)]-\ Since 
the Z-diagram gives no contribution to the nucleon magnetic moments, the 
parameter (3 = 0.645 GeV in Eq. flTUl) can be directly fixed through a fit 
to the experimental values, obtaining /i*^ = 2.87 ± 0.02 (/i^^^ = 2.793) and 
f/^ = —1.85 ± 0.02 (/i^^^ = —1.913). Theoretical uncertainties come from the 
Montecarlo integration of ([7]). 



In Eq. ([HI), the vertex function A(A;i, A;2, /^a)],-;.- ) describes a qqq system 

^ Ion' Son-' 

since < 0, and therefore it cannot be approximated as the one in the 
valence region. It turns out [15] that this NV vertex leads to a contribution 
to the nucleon FF to be interpreted as a transition from \qqq) to \qqqqq) 
Fock components of the final nucleon. In the present calculation, an Ansatz, 



A{ki,k2,k 



2' 1^3) \lk' ) 

^ Ion' Son/ 



, in terms of invariants as the squared free mass. 



M|f(l,2), of the quark pair propagating from the initial nucleon toward the 
final one, and the squared free mass of the system + 3, has been adopted 
(cf diagram (b) in Fig. [1]) 



a5L 



[9i,2n9 



N,3 



|7/2-2 



k 



12 







r 


Pi, 











'111 



where gA,B = {rriA ms)/ [(3'^ + M^{A, B)] and fc^2 = kf+k^. The ratio kt2/PN 
enforces the coUinearity of the spectator-quark pair and the final nucleon, 
while [P'j^ /k'^Y [PN/k^Y controls the end-point behavior of the antiquark- 
leg attached to the nonvalence vertex (with a chosen symmetrical form). The 
powers of gi2 and g^i and the parameter r = 0.17 allow one to obtain a dipole 
asymptotic behavior for the NV contribution. 

In the TL region, where q = P/v + Pr^ ^12 + ^3 = —Pn and ki2 + k'^ = P^, 
after integrating on ki and /c^, one obtains two contributions, with a form 
similar to Eq. ([H]), but corresponding to diagrams (a) and (b) of Fig. [21 In both 
contributions valence and NV nucleon vertexes are present, as a result of a 
transition between the \qqq, qqq) hadronic component of the photon state and 
the NN final state. The nucleon NV vertex is approximated by an Ansatz, 
A^y, analogous to Eq. (ITTi) . built up with the corresponding invariants, e.g.. 
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in the contribution (a) of Fig. [2] one has 



2[^i,2]^[^iV,T2] 



7/2-2 



12 



P 



N 



N 



1.1+ 
L'^ 3 



r r p+nr 

L'^ 3 J 



:i2) 



where the factor 2 counts the possible patterns for gluon emissions. In g^Y2-> 
put in the normahzation factor = 0.500 Gel^[19], and in the denominator 

M2(iV, T2) = (P+ -kt- kt) (P^ - A;i-„ - fc2on) - |ki± + k2xp 

To determine the free parameters Z^, Zy'^j, p and r, a fitting procedure has 
been performed in the SL region, including proton data (/ip G^^/G\j and G^j) 
with Q2 < 10 {GeV/cf and neutron data (G^ and G^^) with < i [GeV/cf. 
We obtained a value x^l datum = 1.7, with a very nice description of the 
data, as shown in Fig. [3l From the fitting procedure we have: i) the ratio 
Z{pj^/Zylj = 1.12, remarkably close to one, and ii) Zb = Zylf = 2.283. In 
correspondence to the previous outcome, the proton charge radius is Vp = 
0.903 ± .004 fm {rf^ = 0.895 ± 0.018 [20]) and -dGl{q^)/dq^ = 0.501 ± 
0.002 (GeV/c)-^ (the exp. value is 0.512 ± 0.013 {GeV/c)~^ [20j). 

The same values for Zb, Zy^ and r (see Eq. ( [T2i) ) are adopted for calculating 
the effective TL form factors, defined as follows, according to experimentalists 
(see, e.g., [21]), 



Glfhq') = ^{\Gf\q^W-v\Gt\<l'W )/(l - V) (13) 

The Z-diagram (higher Fock components) is essential for describing the nu- 
cleon FF, in the adopted reference frame (g^ ^ 0), as in the pion case 
jlOj . In the SL region, it produces the striking feature of a zero around 
~ 9.0 (GeV/c)^ for G^^. Notably, retaining only three sets of data, G^, 
G^ and G^, in the fitting procedure, one gets again the cancellation between 
triangle and pair-production contributions to Gg, and only tiny differences 
from the results shown in Fig. [31 This means that, in our model, the falloff of 
fip G^b/G^m > 1-0 {GeV/cy is enforced by the other three sets of data. 

In the TL region, our calculations are parameter free, and give a fair descrip- 
tion of the proton data, apart the peak at the threshold, which is outside the 
present model due to the absence of the final state interaction. The TL proton 
data clearly show a structure due to the resonances (see Fig. Hj), allowing to 
gather more details on the hadronic content of the photon wave function. In 
particular, the comparison with the most recent data ^1] put in evidence that 
some strength is lacking in our model for ~ 4.5 and ~ 8 {GeV/c)"^ (as for 
the pion [10]). Finally, available TL neutron data are not reproduced by the 
present model, but, even a constant factor of 2 could improve the description 



7 



(cf Fig. m right panel). 

Summarizing, in a frame with ^ Q (that allows a unified treatment of 
both the SL and TL regions), our approach, with only four adjusted param- 
eters {Zb = Zy^, Zyl.j, p and r), is able to describe the nucleon SL FF 
{x^/ datum = 1.7) and to give predictions for the TL ones. A complete analysis 
of the model dependence, as well as a detailed study of the momentum dis- 
tributions of the valence nucleon vertex functions will be presented elsewhere, 
together with a study of nucleon FF in the unphysical region (0 < < 4M^), 
which appears very challenging, but needs a non trivial inclusion of the NN 
interaction [15]. 
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^A^' Pn Pn' (a) Pn Pn' (b) 

Fig. 1. Diagrams contributing to the SL nucleon FF: (a) valence (triangle) contri- 
bution with < fc+ < (i = 1,2,3) and < /c^ + g+ < P'^^; (b) nonvalence, 
pair-production contribution with > fcg" > —q^. A cross on a quark line indi- 
cates a quark on the fe~-shell. Solid circles and solid square represent valence and 
NV vertex functions, respectively; open and shaded circles are bare and dressed 
quark-photon vertexes, respectively. 




Fig. 2. Diagrams contributing to the TL nucleon FF: (a) < + < q^] (b) 
< /c;^ + g+ < Pj^. Symbols as in Fig. [H 
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Fig. 3. Spacelike nucleon form factors vs —q^- Solid lines: full calculation, i.e., sum of 
triangle plus pair-production terms. Dotted lines: triangle contribution only. Data 
from the compilations in [20]. GD{\q^\) = [1 + \q^\/ {Q.7l{GeV / cf)]-'^ . 




6 8 10 
(GeV/c)^ 



q (GeV/c)^ 



Fig. 4. Nucleon effective form factors (see Eq. (1130 for the definition) in the timelike 
region. Solid line: bare + VM. Dotted line: bare term. Left panel: G^^^ ^{q^) / G D^q^) 
vs q^; data from [21]. Right panel: G^^^{q^)/GD{q^) vs q^; data from [5]. Dashed 
line: solid line arbitrarily multiplied by 2. 
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